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Abstract 

We introduce and solve in closed form a simple model of a macroscopic 
body propagating in a one- dimensional resistive medium at temperature T. 
The assumption of completely inelastic collisions between the body and the 
particles composing the medium leads to a resistive force that is opposite 
and proportional to the square of the body's velocity. 
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The topic of macroscopic bodies moving through resistive media, such as 
air or viscous fluids, gives rise to one aspect that students of Introductory 
Physics courses find rather mysterious: The origin of the 'force law' that 
describes the effective force on the moving body as it propagates through 
the resistive medium. The student is usually told that the effective force on 
the body is either proportional to the speed of the body or to the square of 
the body's speed, according to whether the body has, or does not have, a 
small cross-sectional area, or whether it is moving at low or high speeds 1 . 
Under further questioning, the instructor might tell the student that these 
'laws' are based on 'experimental observations' which are difficult to obtain 
analytically. There are basic models, however, that show in a simple manner 
how the energy and momentum exchange between the moving body and the 
particles composing the resistive medium lead to some of these 'force laws'. 
In this article, we present an extremely simplified model that leads to a very 
well-known 'force law': a resistive force that is opposite and proportional to 
the square of the body's velocity: F = —7 V 2 . 

Consider a body, represented by a heavy point 'particle' of mass M and 
initial speed Vq, that propagates inside a one-dimensional medium composed 
of identical point particles of mass m, with m <C M which are in thermal 
equilibrium at temperature T (Fig.l). We will consider here the 'short' time 
scale where the body does not have enough time to reach thermal equilibrium 
with the surrounding medium. The 'brownian motion' case where the body 
is in thermal equilibrium with the medium have been nicely discussed by de 
Grooth 3 . Let us denote by Vj, the velocity of the jth medium particle. Since 
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the medium is one-dimensional the particles can be labelled unambiguously. 
For instance, the particles to the right of the body could be labelled by odd 
values of j, while the ones to the left, by even j values. Because of thermal 
equilibrium the {vj} are random quantities whose values are taken from a 
gaussian distribution of width proportional to the medium temperature T. 
We will assume for simplicity that the body undergoes completely inelastic 
collisions with the medium particles. 

After the first collision we have, because of momentum conservation, 

M V + mvi = (M + m)Vi 

i.e., the speed of the body after its first collision is 

/ M \ / m \ 

Vi= {m^) Vo+ {m^) Vi > 

where v\ denotes the velocity of the medium particle with which the body 
collides first (this particle could come from the left or right of M) . Some time 
afterwards, the body (now with mass M + m) will suffer a second inelastic 
collision from which will emerge with velocity: 

V 2 = ( M ° ) V + ( m ) ( Vl + v 2 ). 
2 \M + 2m) \M + 2m) ' 

where v 2 is velocity of the medium particle who suffers the second collision 
with M , and so on. After n of these collisions, the speed of the body will be 

where we remind the reader that the {vj} are random with (vj) = 0, (vj) = 
kT/m and (...) denotes a thermal average. This implies, 



On the other hand, 



2 M 2 V 2 2mM V / \ m 2 / 

^ n ; (M + nm) 2 (M + nm) 2 \ Y / (M + m) 2 \ ^ ' 
nm/cT 



{Vn? + 

( ///' \ ...... /,-:/ : 



(M + nm) 2 

(V n ) 2 + --—(V n }. (2) 
\ mqVq ) M Vo 

We note that, as the number of collisions tends to infinity (i.e., after a 'long' 
time), (V 2 ) -> (fcT/M V r o)('l4) = (kT/M(n)), where M(n) = M + nm is the 
effective body mass after n collisions. This is nothing else but equipartition: 
M ( n )( v n) ~>kT = m(v 2 ), where v is the velocity of a medium particle. 

If we now assume that p, the density of medium particles per unit length 
is constant, then we can express n asn = p x where x is the distance travelled 
by the body between its first and n-th collision. We are assuming here, as 
in hydrodynamics, that an element of length Ax while 'small' will contain 
a great number of medium particles. By re-expressing n in terms of x in 
(HD, we can express the average velocity of the body after it has travelled a 
distance x as 

(V(x)) = ( ) V 0} (3) 

x v ;/ \M + pm xj ' w 

and the average of the velocity squared as 



The average velocity decreases monotonically with distance. Its explicit time 
dependence can be found from (||): dx/dt = M Vo/(M + pmx), which can 
be integrated to give 

X(t) , 



v ~ v ^ ■ 2(f Ao) (5) 

An 
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where Xq = Mo/(pm) and to = Xq/Vq constitute natural length and time 
scales. Finally, after replacing ([5]) into (|3]), or by direct differentiation of (|5|), 
one obtains 

M> ^1 + 2 (t/t ) 

and 

is the effective body mass as a function of time. In Fig. 2 we show M(t), X(t) 
and V(t), all of which diverge at long times. 

Average resistive force. As the body propagates, it is being hit from front and 
back by medium particles which stick completely to it after colliding. This 
accretion process is rather akin to the opposite process that occurs in the 
propulsion of a rocket engine: instead of expelling matter our body absorbs 
it. One process is the time-reversal of the other. The average effective force 
on the body can be directly computed from (F) = M{t)d(V{t)) / dt. From 
Eqs.(||) and (0) one obtains: 

(F\- MoV ° 1 M 
{F) ~ ~ t l + 2(t/t ) (8) 

which can be recast as 

(F) = - 7 V{tf (9) 

with 7 = pm. 

Another way to compute (F) is to start from conservation of momentum 
during an infinitesimal collision, M (x) V (x)+dM (x)v = (M(x)+dM(x))(V(x) + 
dV(x)). This implies that the instantaneous force on the body is 

„ dV fdM\, Tr N , , 

M s=-[nf) (v - v) (10) 
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where v is random. The average force on the moving body will then be 

<*>--(£) CO ~(£)<v>. 

Since M(x) = Mo + pmx, we now have (F) = —7 (V) 2 , with 7 = pm as 
before. 

Stopping power. The stopping power (5) of a medium is defined by the 
average energy per unit length, lost by a projectile while traversing a resistive 
medium: 

^-{t)=U\ M ^ )2) Y (12) 

From Eqs.(|), (§) and the relations (d/dx)(V(x)) = -(pm/M V ){V(x)) 2 
and M(x) = M V /{V(x)), we have 



,2 



/ feT \ / Mp y 2 
2 ^ M V 2 ) \M + pmx) 
pm (2E - A;T) 2 



(13) 



2M (M V 2 - kT) ' 

Figure 2 shows (5*) as a function of distance traversed inside the medium for 
several temperatures. Note that, since we are assuming the body's initial 
kinetic energy to be higher than the average thermal energy, the body will 
always lose energy to the medium, on average. This energy loss becomes 
smaller and smaller as the body traverses the medium. Only after an infinite 
amount of time, or distance travelled, will the body's average energy loss 
reach zero, where a thermalization process will occur. 

In summary, we have introduced and solved in closed form the dynamics of 
a simple model of a body moving through a resistive medium. We find that 
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the effective resistive force is opposite and proportional to the square of the 
body's speed 2 . 
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Figure Captions 



FIG 1: Macroscopic body of mass M propagating inside a one-dimensional 
resistive medium composed by identical particles of mass m « M in ther- 
mal equilibrium at temperature T. 

FIG 2: Effective body mass, average velocity and distance travelled as 
a function of time, for body moving through our resistive medium (So = 
P mVi/2). 

FIG 3: Stopping power of the one-dimensional resistive medium as a function 
of the distance traversed by the body, for several medium temperatures. 
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FIG. 3 



